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We propose a kinetically driven mechanism based on the breaking of the spatial reflection sym¬ 
metry to describe the magnetic moment and the torque observed by Lu-Li et al. (Ref. [J) for the 
LaAlOs/SrTiOa system. We find that the itinerant electrons are in a zero helicity state and predict 
the existence of charge inhomogeneities that cross the interface at constant rate. There is mass and 
thickness anisotropies between the two sides of the interface. 


PACS numbers: 74.20.De,74.25.Ha,75.70.Kw,75.70.Tj 

Introduction. - The discoverj*2i^ of a high-mobility 
electron gas at the interface between the bulk insula¬ 
tors non-magnetic LaAlOa (LAO) and SrTiOa (STO) 
unveiled unexpected properties^, such as the coexistence 
of magnetism and superconductivity. The Rashba spin- 
orbit interaction^ plays a pivotal role^ acting on both the 
superconducting and the normal-state transport proper- 
tiesi. There are signals of an underlying ferromagnetic 
order— that appears as magnetic dipoles, but with no net 
magnetization^. Torque magnetometryi shows the exis¬ 
tence of a magnetic moment parallel to the interface for a 
nearly perpendicular applied magnetic field. To reconcile 
the observed perpendicular (axis 3) and parallel mag¬ 
netic moments, theoretical proposals of elaborate mag¬ 
netic patterns have been made, based on spiralsiS, and 
skyrmions^, obtained from the ad hoc assumption of the 
Rashba term. A general two dimensional (2D) free en¬ 
ergy has been suggested to describe the long-wavelength 
magnetisroi^ to find a lattice of skyrmions similar to he- 
limagnets. In this letter we propose that the magnetic 
moment and the torque observed in the LAO/STO sys¬ 
tem are kinetically driven, which means that they stem 
from the standard 3D kinetic energy of itinerant electrons 
with no need for extra terms. Nevertheless the presence 
of an interface breaks the reflection symmetry and gives 
rise to the Rashba term, which is then found to be con¬ 
tained in the 3D kinetic energy. As a consequence of 
this breaking the electronic confinement to the interface 
is quasi two-dimensional, the itinerant electrons are in a 
zero helicity state and there currents crossing and enter¬ 
ing the interface constantly. The existence of a lattice of 
vortices and skyrmions in the interface follows from this 
zero helicity condition (ZHC). To describe the magnetic 
moment and the torque data of Ref. |l| there must exist 
anisotropyi^, between the LAO and STO sides of the in¬ 
terface. The mass of surface carriers has been detected 
since long ago in superconducting thin films and wires 
deposited on a SrTiO^^. 

The torque in the high-Tc superconductors is known 
to be kinetically driven and is well described by the phe¬ 
nomenological London anisotropic theoryi^. The Lon¬ 


don kinetic energy captures well the energetic unbalance 
caused by a supercurrent circulation that preferably re¬ 
mains along the easy mass plane, namely, the supercon¬ 
ducting layers, for a tilted external field. As a result there 
is magnetic moment not oriented along the applied field 
and a consequent torque whose measurement can unveil 
new properties of the vortex state. The temperature de¬ 
pendence of the anisotropy has been used to give evi¬ 
dence of two band superconductivity in the layered com- 
poundsi^. Torque measurements helped Lu-Li et to 
claim that there is a rigid London order parameter above 
Tc despite the loss of the Meissner effect in the high-Tc 
compound Bi2Sr2CaCu20s- Theoretically it has been 
shown that the torque can show signals of vortex-vortex 
attraction, known to occur in the low tilted field regime of 
the high-T(i^. The unbalance caused by the coexistence 
of vortices of different lengths in granular superconduc¬ 
tors leads to a torquei^. In this letter we claim that the 
torque observed in the LAO/STO system is also kinet¬ 
ically driven. Fitting the magnetic moment and torque 
data of Ref. [l| to the present theory gives information 
of the mass anisotropy and the effective thickness of the 
LAO and STO sides of the interfacei^. 

The theory. - Interestingly, according to Ref. [l|, the 
magnetic ordering and the superconducting state coexist 
and the torque remains nearly the same far above the su¬ 
perconducting critical temperature. We take this as an 
evidence of a kinetically driven torque, namely, that nor 
the condensate energy, below Tc, neither the interaction 
energy, above Tc, are relevant to the torque. Therefore 
our starting point is the simplest possible order param¬ 
eter theory composed by the sum of the 3D kinetic and 
field energies, F = Fk -\- Ff, where Ff = takes 

the contribution of the local field h generated by the cir¬ 
culating currents. Notice that this starting point is ex¬ 
actly the same one taken to describe the torque of the 
high-Tc superconductors^^. However a fundamental dif¬ 
ference assumed here is the presence of two order param¬ 
eters, each one associated to a distinct spin component. 
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FIG. 1. (color online) Measured magnetic moment and torque 
data fitted by the present theory 
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(b)The local field (5/i||(0'*') 
(parallel and just above the 
interface). 


(c)The local field (5/i||(0“) 
(parallel and just below the 
interface.) 


FIG. 2. The local magnetic field without the applied one is 
shown in the unit cell (xi/Z/,X 2 /T), 


The 3D kinetic energy is well known to be given by, 


where 'F is the two-component order parameter. The 
minimal coupling is given by D'f. = ymo/rn^Dk, 
D = (?i/i)V — {q/c)A, where m^k) are the mass pa¬ 
rameters along the major axes, parallel to the inter¬ 
face (fc = 1,2) and perpendicular to it (k = 3). Uni¬ 
axial mass symmetry is assumed, such that the paral¬ 
lel masses are equal in both the LAO and STO sides, 
TO(i) = 771(2) = ” 7 , but the perpendicular masses are dis¬ 
tinct in each side, and given by 777 .( 3 ) = ^ ^ind 777 ( 3 ) = -V, 
respectively. The thickness of the LAO and STO sides 
are d and d, respectively, such that the active volume is 
5V = A{d + (£), A being the area of the interface. There¬ 
fore the integration volume comprises the LAO and STO 
volumes separately, (•••) = (•• •)x 3 >o + (• • •)a: 3 <o, where 
{■■■)x 3 >f) = dx^ {■■■)/{Ad) and {■■■)x 3 <o = 

f-d^^3 {■■■)/{Ad). 

Remarkably the 3D kinetic energy admits the follow¬ 
ing distinct but equivalent formulation^^, suitable for the 
breaking of the spatial reversal symmetry. 


Fk = 


27770 


2 


^rt 


CT- 


(ct X D') 


-I- -h • 
c 


^ + c.c. 


( 2 ) 


where V' = y^rn^Jm(i)Vi and ct' = m^/ 

7770 being an arbitrary parameter. The Pauli matrices 
are Gi and i ^ j ^ k take values among 1 , 2 , and 
3. Eq.([5]) splits the 3D kinetic energy into three con¬ 
tributions. The last term, proportional to ct x D', is 
the Rashba term, only present at the edges, which here 
corresponds to just {xs = 0 +) above and {xs = 0 “) 
below the {x^ = 0) interface. Along the interface we 
assume periodicity defined by a unit cell. The helicity, 
u-D', is a pseudo scalar, and consequently, the imposition 
that IF makes it vanish triggers the breaking of the spa¬ 
tial reversal symmetry. Interestingly the ZHC yields the 
minimum electromagnetic energy since it exactly solves 
Ampere’s law, and fully determines the local magnetic 
field h = V X A. To see this consider the electromag¬ 


netic current, known to be J = {q/mo) 


+ c.c. 


in the traditional formulation of Eq.([T]). In our al¬ 
ternative formulation of Eq.Q, the current is J = 


{q/2mo) I (ct.ZI'If) 
one obtains that, 


-I- c.c — hV‘ 


' X 


Thus 


G ■ D' = 0, and (3) 

5h = h-H,,, = -y^-^^G''l>, (4) 

7770 C 

where 5h is the local field without the external applied 
field. A most important remark is that h, obtained from 
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(a)—V|| • Js (the charge rate density) 
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(b) Js (the superficial current) 
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(c)l’Ifp (the order parameter density ) 


FIG. 3. Properties at the interface = 0. 


Eq.Q, is a solution of Maxwells equations, which means 
that V • h = 0 . Thus h stream lines pierce the inter¬ 
face twice or none. The first case implies in the exis¬ 
tence of a superficial ( 2 D) current density Jg at the in¬ 
terface since ^3 x [h( 0 +) — h( 0 “)] = 47rJs(0)/c. This 
current gives rise to the charge inhomogeneities along 
the interface since Vy • Jg ^ 0 . We obtain the mag¬ 
netic moment and the torque under the two assumptions 
of a (i) small order parameter T Ri 0(e) and of a (ii) 
small tilt angle 9 between the applied field and the inter¬ 
face. The parameter e controls the smallness of the fields, 
for instance, from Eq.Q it follows that 5 h Ri O(e^). 
The perpendicular and parallel components of the ap¬ 
plied field with respect to the interface, are given by 
H 3 = Hext cos 9 and iJ|| = Hextsm9, respectively, such 
that Hext = H3X3 -|- where ey is a vector parallel 

to the interface. We solve Eqs. m and (HJ recursively to 
order O(e^), which means that 'h is obtained from Eq.Q 
for H3X3, with iJy discarded because of the assumption 
(ii). Next the inhomogeneous field 6 h is obtained from 
Eq.Q using the knwon Further iterations of Eqs. ([ 3 ]) 
and (HI) are not necessary to this lowest order O(e^). 

The torque directly follows from t = m x H^xt, once 
known the magnetic moment rh = VM, where 

is the magnetization. The magnetic induction is B = (h), 
and comparison of Eq. (| 3 ]) with B — Hext + 47rM yields 
the above expression. Another way to obtain the torque 
is from the free energy, r = —VdF/d 9 , and as we are 
restricted to the lowest order O(e^), only the following 
terms can contribute to it. 

^+ ( 6 ) 

Next we show that the same torque can be obtained ei¬ 
ther from this free energy or the magnetic moment di¬ 
rectly derived from Eqs. and (| 3 |). It can be shown 


that the last term of Eq.ljS]) is proportional to H^, and 
so, it does not contribute to the torque under assump¬ 
tion (ii). A way to see this proportionality is to con¬ 
sider that this last term vanishes for a spatially homo¬ 
geneous state, namely Tl'k constant, and this is only 
achieved for H3 = 0 . Therefore the free energy becomes 
F r; —Hext ■ M that can be reduced to F r: —H\\ ■ My, 
where My = —{qh/2mQc){^'^a'^^"^), CTy = cr[xi + a{xi. 
There is no net magnetization perpendicular to the inter¬ 
face, M3 = 0, in the present approach, which is in^ree- 
ment with the experimental observations of Ref. As 
shown below My only depends on H3, and so the remain 
non-vanishing contribution under assumption (ii) comes 
from t/V Ki —dHji/d9 ■ Mji = i/3(g7i/2moc)ey • ('l>l(j|| 4 '), 
thus rendering the same torque of Eq.(l 5 ]). 

The magnetic moment parallel to the interface. 
- The order parameter that satisfies the ZHC is obtained 
under the choice of gauge A3 = 0 in Eq.(l 3 l), and is. 






'tpn(xi,X2) \ 
lftjV'n-i(a;i,X2) J ’ 


( 7 ) 


It contains the Landau level functions, ■i/'n, that ap¬ 
pears twice in this series with the exception of ipo ■ They 
are normalized to /^gjj(d^a;/L^)|'(/i„p = 1, where the cell 
length is L = The ZHC does not fully deter¬ 

mines the order parameter since any set of coefficients c„ 
provide a new solution of Eq.(I 3 ]). We expect that resid¬ 
ual interactions and higher order terms in the free energy 
determine the coefficients c„ but we do not carry this pro¬ 
cedure here. We take that the torque is not significantly 
sensitive to the free energy minimization, similarly to the 
case of the high-Tc superconductor case. There the ex¬ 
perimental torque cannot determine the optimal param¬ 
eters of the vortex lattice. 

The order parameter decays away from the interface 
and is distinct in each side of it. This decay is determined 
by a = r]{m/MY^‘^ {x^ > 0 , LAO) and a = r]{m/, 
{x3 < 0 , STO), respectively where ry = ^/WqJHcHZ- 
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Therefore the magnetic field controls both the periodicity 
at the interface, and also the evanescence perpendicular 
to it, through L and ij, respectively. Indeed the exponen¬ 
tial decay described by Eq. dZD sets a quasi 2D behavior 
away from the interface, that can be either 2D or 3D, 
and is controlled by a, and so, by rj and the anisotropy. 
The magnetization is given by. 



system, total area of the interface is A r; O.dlO^^nm^, 
and the height is d + d = 4.3 nm. We notice that the 
maximum magnetic moment seen in Fig. (HD corresponds 
to the parameters r and r of order of a few units, thus 
showing that the exponential behavior of Eq. dzD sets the 
transition form 2D to 3D behavior. From this theory 
it follows that a qualitative value for the local field is 
(5/i|| Ri 47rm|| for a given H^xt- For instance at the lowest 
applied field of H^xt = 5 mT the measured moment^ is 
of TO|| = 4.35 and this gives that J/iy r; 70 mT. 


where /is = Tiq/2mc, and we find that the two terms 
contribute oppositely, 

OO 

= -ix\\ ^ c„*c„+i/(„) (r) - c.c. (9) 

n—1 

OO 

(TV||4')^^<0 = ixi\ Cn*Cn+lf(n) {r) - C.C. (10) 

n —1 

where aiy = (ii -|- ix 2 ) and only the ratios r = 
{MImY^'^dlr] and f = {M/mY/’^d/rj matter for such 
averages and enter through the function /(„) (z) = 
[1 - exp (-Zn)]/z„, Zn = {^/nTT + y/n) z. No¬ 
tice the general property, a1 || d, d] = 

—M|| [(■^)^^^, d, (^)^^^,d] which promptly shows that 
M||(—ids) = —M||(id 3 ) since ^ —ids is equivalent 
to reverting the LAO and STO sides, namely, d d 
and M M. Thus the present approach is in agree¬ 
ment with the reflection symmetry property observed 
in Fig.([T|). The asymptotic limits of My ^ and 

M|| ^ l/y/Hs are a property of Eq.([8D for small and 
large fields, respectively and if added to the fact that 
M|| does not changes sign for ids > 0 , shows that the 
magnetization must reach a maximum, as seen in Fig. (HD. 

Properties and comparison with the measured 
torque. - The above general properties for the mag¬ 
netic moment and torque are valid for any set of coef¬ 
ficients Cn in Eq. HD and allow for fitting of the data of 
Ref. [H. As pointed before ipn appears twice, in consecu¬ 
tive doublets, therefore the simplest set able to capture 
the essence of this expansion corresponds to ci = C 2 = e 
real and c„ = 0 for n < 3. With respect to the choice 
of fitting parameters, the thickness of the LAO layer is 
^ 5 unit cellsii^i^ (5a m 2.0 nm, a m 0.4 nm). Un¬ 
der these conditions we find for the best fitting the pa¬ 
rameters d = 1.8 nm, d = 2.5 nm, {MjmY^'^ = 8.0 
and {MjmY^'^ = 8.5 and is shown in Fig. HD- From 
it we obtain that e « 9.4nm“^/^. From fitting it fol¬ 
lows that Ue^ = 1.610^® and V fv l.SlO^^nm^. The 
volume is obtained by assuming 0.3 — O.d/r^ moments 
per LAO/STO cell area a^, according to Ref. |l|, and 
since there are approximately 10^^ fj,B moments in the 


Figs.dD) and ([3]) show properties of the state at Hgxt = 
0.5 T. The local magnetic field near to the interface is 
shown in Fig. HD- Fig.( |2(a)D reveals positive and nega¬ 
tive puddles of < 5/13 at the interface, normalized to arbi¬ 
trary units. This shows the existence of closed 6h stream¬ 
lines encircling the interface, confirmed by Figs. (2(b)) 
and (2(c)) which show dh\\ immediately above and be¬ 
low the interface, respecitively. The point (0.5, 0.8) is 
the skyrmion core as Sh stream lines cross from one 
side to the other of the interface to return elsewhere 
in the unit cell. There is topologically stability for the 
skyrmion as the number, Q — (l/47r) \[dh/dxi) x 

{dh/dx2)\ - h d'^ X, is an integer {Q = —2). There is 
also vorticity carried in the phases of The 3D and 
2D currents, J and Jg, respectively, render a total di¬ 
vergenceless current at the interface too, which means 
that V|| - Jg + </ 3 ( 0 '^) — J 3 { 0 ~) = 0 and is interpreted 
as V|| - Js + da/dt = 0. Charge crosses the interface 
at constant rate da/dt defined by the charge conserva¬ 
tion forming positive and negative puddles as seen in 
Fig.(3(a)). We obtain from the present model that qual¬ 
itatively da/dt ~ 5h\\c/L. Under the assumption that a 
charge unit e crosses a nm^ area this expression defines 
a rate that falls in the THz range. Fig.(3(b)) depicts Jg 
and shows a circulation around the skyrmion center at 
the point (0.5,0.8). However this does not corresponds 
to a zero of the order parameter density, as shown in 
Fig.(3(c)). This is because the order parameter has two 
components and only one of them vanishes at this point 
meaning that indeed there is a vortex there, but asso¬ 
ciated to only one of the two available phases. In con¬ 
clusion we have shown here that the ZHC explains the 
magnetic moment and torque data of Ref. [l| from a ki- 
netically driven mechanism. 
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